Anomalous temperature dependence 
of the order parameter of a 
superconductor with weakly 
correlated impurities 



It is shown that weak correlations between pair-breaking impuri- 
ties in superconductors influence the temperature dependence of the 
order parameter within the Ginzburg and Landau region if correla- 
tion radius of impurities R is greater than the coherence length of 
the superconductor £o . Dependence of a square of the average order 
parameter on the temperature difference T c — T changes its slope in 
a region £o \JT C /{T C — T) ~ R . Influence of correlations of impurities 
on other thermodynamic properties of superconductors is discussed. 

Impurities make the condensate of Cooper pairs spatially nonuniform. 
Manifestation of this non-uniformity is particularly strong in unconventional 
superconductors and in conventional superconductors with magnetic impu- 
rities. In a vicinity of the transition temperature T c free energy of such 
superconductor can be written as Ginzburg and Landau functional with 
coefficients, which are random functions of coordinate. For a scalar order 
parameter \&(r) 
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Abstract 




According to the analysis of Larkin and Ovchinnikov [TJ the most strong 
effect on the average order parameter have fluctuations of the coefficient 
a(r) i.e. fluctuations of the local transition temperature T c (r) : a(r) = 
a(T — T c (r)) . Because of these fluctuations the temperature dependence 
of the average order parameter ( 1 J r ) deviates from the linear dependence 
(\]>) 2 r*j (T c —T)/T c characteristic of the uniform superconductor and becomes 
singular in a narrow interval below the superfluid transition temperature 
(T c — T)/T c ~ (^ /0 4 (V( n ^o) 2 ) , where I is the mean free path and n - 
the density of impurities. The linear dependence is preserved outside of 
this region. This result is obtained under assumption that impurities are 
not correlated. It has been shown recently, that correlations with a radius R 
which is greater than £ strongly affect suppression of T c by impurities 013]. 
Such situation is realized for the superfluid 3 He in aerogel. Experimental 
data for thermodynamic properties of superfluid 3 He in aerogel, such as the 
square of the Leggett frequency Q 2 L and p s [U E] deviate from the linear 
dependence on T c — T in a much wider interval of temperatures than that, 
estimated on a basis of Ref. [Tj. The observed dependence of these quantities 
bends upward when T c — T increases. The aim of this paper is to show 
that this anomaly can be interpreted as the effect of correlations among the 
impurities. The effect of correlations is general in a sense, that presence 
and a character of deviations do not depend on a particular form of the 
order parameter. For demonstration of the effect the simplest example of a 
superconductor with the scalar order parameter ^(r) will be considered. 

Let us follow the argument of Ref. [1] with the modifications required by 
the presence of correlations. Free energy (1) can be rewritten in terms of the 
dimensionless order parameter ip{r) = ^(r)/^o, where \l/o is the absolute 
value of the order parameter for pure superconductor at T = obtained by 
extrapolation of linear dependence of |\I/| 2 on (T c0 — T) from the transition 
temperature of the pure superconductor T cQ : 

F s {^(v)} = F n + T c0 Ac J {[r-^(r)]|^(r)| 2 + ^(r)| 4 + C 2 |V^(r)| 2 }rf 3 r. (2) 

Other notations here are: Acq - the specific heat jump in the pure super- 
conductor, r = (T — T cl )/T c0 , T cl = (T c (r)) , so that the relative fluctuation 
of local transition temperature rj(r) = (T c (r) — T cl )/T c0 , it vanishes after 
averaging. The elasticity coefficient £ 2 = -^rp^o m the BCS theory. In these 
notations Ginzburg and Landau equation reads as: 

[ T -77(r)]^ + V#| 2 -£ s 2 AV> = (3) 

At T < T c the long-range order is established, characterized by the average 
order parameter (ip) . The angular brackets denote averaging over ensemble 
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X » X + X ' X- X + X ' X ' X ' X + ... 



Figure 1: 

of impurities. Solution of Eq. (3) can be sought in a form ip(r) = + 
x( r )) • When r](r) is small x( r ) is small too, except for the mentioned 
above temperature region close to T c , where (x(r)x(r)) diverges. Averaging 
Eq.(3) over ensemble of impurities renders expression for (tp) 2 in terms of 
the average products (rjx) and (xx) '■ 

W 2 = {vx) 7\ (4) 

The fluctuating part x( r ) can be found from the linearized equation (2): 

(r - 77(r) + 3(^) 2 ) X - t^X = ~ (vx)- (5) 
Its solution can be formally written in terms of the Green function G(r, r') : 

X (r) = / ^,0(77(0 -(vx))d 3 r'. (6) 

Re-writing this solution in the momentum representation we can find the 
averages, entering Eq. (4): 



d 3 k 



(V(r)x(r)) = J (v(-VG(kX) V (k>)) — — [l + / (,(-k)G(k, 0)) ^ 

(7) 

The average in the numerator can be found by term-by-term averaging of 
the diagrammatic series Fig.l. 

To every arrow here corresponds the unperturbed Green function 
G (k,k') = (27r) 3 5(k - k')[r + 3(^) 2 + gk 2 ]" 1 and to a cross 77(1*1 - k 2 ) , 
ki is the in-coming and k 2 - the out-going momenta. Comparison of the 
averaged diagrammatic series Fig.l with that for the average Green function 
relates (rjGrj) to (G) : 

(^) = G 1 [(G)G 1 -1], (8) 
which in its turn can be expressed in terms of the self-energy E(k, r) : 

(G(k)) = [r + 3W 2 -S(k,r)+C 2 A: 2 ]- 1 (9) 
Using analogous argument for evaluation of (r]G) we find eventually that 

( X 77>=E(0,t). (10) 
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In a principal order on the perturbation r](r) 

l ' J 7 r + 3(^) 2 -S(k 1 ,r)+e^(2vr) 3 ' 1 J 

Of practical interest is a situation when r/(r) = X) ^ ( r— r o) j where ^'(r— 
r n ) is a contribution of individual impurity situated at the position r a . 
In that case the correlation function (^(-k)r](k)) = n|?7( 1 ^(k)| 2 5'(k) , where 
= (J2b exp[ik(r b — r )]) is the structure factor. For non-correlated 
impurities only one term with = r a contributes to the sum. Effect of 
correlations is contained in SQs) = S(k) — 1 , which is Fourier transform 
of the correlation function in the coordinate space v(r) . When impurities 
are correlated on a distance ~ R the S(k) is enhanced for k ~ 1/R. For 
estimation of the effect of correlations a model expression (" 8 -model" [3]): 
S(k) = 2n 2 R 2 nv(0)5(k —1/R) can be used. In particular, one can show 
that corrections to the principal expression for £(0,r), given by Eq.(ll) 
contain extra factor R 2 /^qI = e. In what follows we assume, that £ C 1 
and treat effect of correlation as a perturbation. In a principal order on e in 
the denominator of the integrand in Eq. (11) we can set (ip) 2 = S(0, r) — r 
and take S(k, r) « S(0, r) . Then Eq. (11) determines a function S(0, r) . 
In this equation a combination u = E(0, r) — r is taken as a new variable. 
In the expression (11) for S(0, r) we separate part which is finite at u = . 
As was discussed before [2j [3], it determines the second order correction to 
the shift of T c . The new transition temperature T c2 is then T c2 = T& + 

T c o / ^kgfcl jj^js ■ It is convenient to introduce a variable t , which 
counts temperature from T c2 : t = (T — T c2 )/T c0 . In these notations Eq.(ll) 
renders relation between u and t : 

|r/( 1 )(k 1 )| 2 S(k 1 ) d 3 h 



u 



1 + 2 



n 



-t. 12 



e s kl[2u + Qk 2 ] (2tt) 

In terms of w (^) 2 = tt/(l+3(xx)) . When 5(ki) is strongly peaked at k\ ~ 
1/i? the dependence of u on t changes its slope in a region u ~ {^ S /R) 2 , 
i.e. when f(T) ~ i? and ^(T) is defined as Wv 7 ^- At u » Hs/Rf 
asymptotically u ~ (T cl — T)/T c0 , i.e. it depends linearly on a distance from 
the average transition temperature T cl . In the opposite limit u (t; s /R) 2 
u is linear in t = (T c2 — T)/T c0 , with a different slope. The relative change 
of the slope is on the order of e 2 . Using for S(k) the 5 -model we have: 

U= l + 2n 2 \r]W(0)\ 2 v(0)(R/^' ^ 

If impurities, like aerogel, have a tendency to form clusters, v (0) > and the 
slope of u(t) at ?j < (t; s /R) 2 is smaller than at u ^> (£ s /i?) 2 , so the u(t) 
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bends upward. Together with u(t) changes its slope {^) 2 (t) = u/(l+3(xx)) ■ 
The average (xx) i n the denominator is on the order of e 2 : 

r n[^)(k)pg(k) _gfc_ 
UX/ 7 [2« + £ s 2 A; 2 ] 2 (2tt) 3 ' 1 J 

This correction does not influence asymptotic of the dependence of {ijj) 2 
on t at u ^> (£, S /R) 2 , but at w (t; s /R) 2 it increases the change of the 
slope. Due to this correction physical quantities, which depend on averages 
of different powers of the order parameter will have different changes of the 
slope. E.g. the NMR frequencies are proportional to (ip 2 ) = (0) 2 (l+(xx)) — 
u/(l + 2 fox)) . For the model S ~ S(k - l/R) at «« (Zs/R) 2 ■ 

^ = l + An 2 \r]m(0)\ 2 v(0)(R/Q^ ^ 

Further thermodynamic properties can be found as derivatives of the av- 
erage free energy over T . Using the expression F s — F n = — Tc0 ^ c ° (J \^\ 4 d 3 r) 
and relation (l^l 4 ) = (V ; ) 4 (l + ®{xx)) we arrive at: 

j? rp T c0 AC 2 

F s - F n = u , (16) 

i.e. a gain of the free energy is proportional to u 2 and not to t 2 . The specific 
heat jump, following from Eq. (16) is: 

AC, = AC T c0 T c2 (^) . (17) 

For the 5 -model AQ = AC g(l + 2n 2 \r]W(0)\ 2 v(0)(R/Z s ) 4 )- 2 . There is 
extra suppression of the jump due to correlations. 

So, we conclude, that the anomaly of temperature dependence of ther- 
modynamic properties of superfluid 3 He in aerogel is qualitatively the same 
as that found in the considered example of a superconductor with a one- 
component order parameter and correlated impurities. The quantitative 
comparison of the obtained results with the data for 3 He would not have 
sense, because superfluid 3 He has a multi component order parameter. This 
opens a possibility for correlated impurities to interact with several different 
modes of fluctuations of the order parameter. For a quantitative description 
of the anomaly all of these modes have to be taken into account. Neverthe- 
less a qualitative estimation of the correlation radius based on the relation 
£(T) f» R at a temperature of a change of the slope is close to other estima- 
tions. 
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Because of its universal character the anomaly can occur in superconduc- 
tors with different order parameters and can be considered as an indication 
that impurities are correlated and correlation radius is greater than £o • A 
distance from T c at which dependence of (ip) 2 on T c — T changes its slope 
renders estimation of the correlation radius. 
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